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Niels Bohr in the early stage of his career developed a nonlinear theory of fluidic surface oscillation 
in order to study surface tension of liquids. His theory includes the nonlinear interaction between 
multipolar surface oscillation modes, surpassing the linear theory of Rayleigh and Lamb. It predicts 
a specific normalized magnitude of 0.41677^ for an octapolar component, nonlinearly induced by a 
quadrupolar one with a magnitude of r) much less than unity. No experimental confirmation on 
this prediction has been reported. Nonetheless, accurate determination of multipolar components is 
important as in optical fiber spinning, film blowing and recently in optofluidic microcavities for ray 
and wave chaos studies and photonics applications. Here, we report experimental verification of his 
theory. By using optical forward diffraction, we measured the cross-sectional boundary profiles at 
extreme positions of a surface-oscillating liquid column ejected from a deformed microscopic orifice. 
We obtained a coefficient of 0.42 ±0.08 consistently under various experimental conditions. We also 
measured the resonance mode spectrum of a two-dimensional cavity formed by the cross-sectional 
segment of the liquid jet. The observed spectra agree well with wave calculations assuming a 
coefficient of 0.415± 0.010. Our measurements establish the first experimental observation of Bohr’s 
hydrodynamic theory. 

PACS numbers: 05.45.-a,42.25.Fx,47.85.Dh 


I. INTRODUCTION 

Optofluidics deals with a synthetic system of optical 
and fluidic elements BHa, utilizing advantages of both ar¬ 
eas [3]. As one of optofluidic examples, deformed micro¬ 
jet cavities are two-dimensional optical resonators sup¬ 
porting whispering-gallery-mode(WGM)-like cavity res¬ 
onances [ 3 . A deformed microjet cavity is formed by a 
cross-sectional segment of a fluidic microjet, which un¬ 
dergoes surface oscillation as shown in Fig. [l] It pro¬ 
vides not only a high quality factor owing to its clear 
and smooth surface but also high output directionality 
based on internal ray and wave dynamics. These fea¬ 
tures can be tuned in real time by controlling the flow 
rate or ejection pressure. Deformed micro jet cavities are 
known to be a versatile platform for studying quantum 
chaos, prototyping highly efficient light sources as well as 
various photonic devices EHin]- 

As is often the case with tunable optofluidic compo¬ 
nents, the optical characteristics of a deformed micro¬ 
jet cavity is sensitive to the detailed shape of its sur¬ 
face boundary. Even for non-optofluidic application such 
as optical fiber spinning and film blowing m , accurate 
determination of multipolar components is required for 
consistent material processing. For optical cavity ap¬ 
plications, in particular, it is necessary to measure the 
boundary profile of this cavity as accurately as possible. 
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However, it is generally not easy to measure the accurate 
shape of a fluidic object. Liquid surface exhibits specular 
reflection and thus make the conventional optical trian¬ 
gulation technique inapplicable [12| . Moreover, the con¬ 
tinuous columnar shape of the microjet makes top-view 
imaging of its boundary simply impossible. 

It may seem that these difficulties in measurement 
might be easily compensated for by a theoretical ap¬ 
proach. It is particularly because there exists a well- 
established hydrodynamic theory to describe the surface 
oscillation of a liquid jet, mostly owing to Lord Rayleigh 
and Sir Lamb [HHH]. However, this theory, based on 
the first-order approximation, cannot explain some de¬ 
tailed experimental facts from spectroscopic observations 
as to be discussed below. It was Niels Bohr at his early 
career who first extended the surface oscillation theory 
by including nonlinear interactions m- The extended 
theory and accompanying surface tension measurement 
of his own based on the lowest-order surface oscillation 
period enabled him to win a gold medal |I6j in an aca¬ 
demic contest of the royal Danish academy. However, the 
surface oscillation period does not reveal the multipolar 
surface-oscillation amplitudes, which is the hallmark of 
Bohr’s theory. So far, direct experimental verification 
of the multipolar surface oscillation amplitudes has not 
been reported. 

In this paper, we present two independent experimen¬ 
tal studies to determine the multipolar surface-oscillation 
amplitudes in the boundary profile of a deformed mi¬ 
crojet cavity. One is a non-destructive surface recon¬ 
struction experiment based on optical forward diffrac- 
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tion. The other is the spectroscopy of the cavity reso¬ 
nances compared with wave calculation results for various 
cavity boundary profiles. We found the results of both 
studies exhibit good agreement with Bohr’s prediction 
in nanometer-level accuracy. Our work thus marks the 
first experimental verification of Bohr’s hydrodynamic 
theory. Moreover, precise knowledge of our microjet cav¬ 
ity boundary allows us to predict optical and dynami¬ 
cal properties of cavity resonances in advance. Quality 
factors, output directionality and intermode interactions 
can be predicted beforehand to guide actual experiments 
for various photonics applications. 

This paper is organized as follows. In Sec. [llj with the 
linear theory of Rayleigh and Lamb applied to our liq¬ 
uid jet, the basic properties of the jet are discussed. We 
then explain how optical cavity resonances are observed 
for our jet and discuss their inconsistency with the cross 
sectional shape of the jet predicted by the linear theory. 
The inconsistency is attributed to small shape perturba¬ 
tions, based on high sensitivity of intracavity ray dynam¬ 
ics and consequent cavity resonances on such small shape 
perturbations. In Sec. m the principle of our optical for¬ 
ward diffraction technique and the experimental result on 
microjet surface profile measurement are presented. In 
addition, the wave calculation results assuming the sur¬ 
face profile obtained from the forward diffraction experi¬ 
ment are presented and shown to be consistent with the 
observed cavity resonances. In Sec. IV Bohr’s nonlinear 
theory, its limitations and our extension are discussed. 
We then show Bohr’s theory is adequate for our liquid 
jet under the current experimental condition. In Sec. 
we summarize our results with remarks on possible ap¬ 
plications and future prospects. 


II. APPLICATION OF LINEAR ANALYSIS TO 
A LIQUID JET 



FIG. 1: Examples of the fluidic surface oscillation due to sur¬ 
face tension, (a) Coffee poured from a pot. (b) Water from a 
faucet, (c) A 3-D model of a vertically ejected liquid jet. The 
pattern is stationary while the liquid moves in 2 ; direction. 
Extreme positions of surface oscillation with stationary cross 
sections are indicated as Dl, D2, D3, and so on. 


A. A liquid jet to be analyzed 

A liquid stream from a non-circular orifice shows semi- 
periodic surface oscillations due to surface tension as a 
restoring force. Such oscillations are easily observed in 
our daily lives as illustrated in Fig. The deformed 
boundary shape of our microjet cavity is also due to this 
surface oscillation. In our setting, a liquid jet is verti¬ 
cally ejected from a near-elliptical orifice [S] as shown in 
Fig. [^c). When the ejection pressure is held fixed, the 
surface shape of the jet becomes stationary. The major 
and minor axes of the horizontal cross-sections are re¬ 
peatedly exchanged as we move along the vertical axis. 
This phenomenon is known as ‘axis switching ’[17]. More¬ 
over, the amplitude of surface oscillation decays and thus 
the cross-sectional shape would gradually converge to a 
circle. 

A partial side view of the jet is shown in Fig. [^a) and 
the microscope image of the orifice is shown in Fig. ^b). 
The shape of the orifice is decomposed into quadrupolar 


(labeled by 2) and octapolar (labeled by 4) components 
mostly. The decomposition serves as one of the initial 
conditions in the analysis below. From the microscope 
images of the jet like Fig. [^a), we can easily identify 
the directions of major and minor axes of cross-sectional 
shape. Moreover, the axis-switching is also clearly ob¬ 
served. The average period of the axis switching for the 
jet is 270 ± 10 /rm in Fig. [^a). 

B. Linear Analysis 

According to the linear theory of Rayleigh and Lamb 
[I^15j . the mth mode of surface oscillation can be ex¬ 
pressed in the cylindrical coordinates as 

rm = arjm cos{m9) cos(kmZ + ^m)e~''^^"' (m = 2, 3,...) 

( 1 ) 

where a is the mean radius, 77 ^ is the relative amplitude, 
km = I is the wave vector in 2 direction with A^ 
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FIG. 2: Actual shape of our liquid microjet. (a) Microscope 
image of the liquid jet column seen from the side. Mean ra¬ 
dius is 15 fj.m. The bright region of about 3 fim thickness 
is the segment excited by a pump laser beam to form a two- 
dimensional microcavity, (b) Microscope image of the orifice 
filled with liquid, (c) Multipolar decomposition of the bound¬ 
ary profile of the orifice of (b). The origin is chosen so as to 
make the dipolar component vanish. In this decomposition, 
quadrupolar (denoted by 2) and octapolar (denoted by 4) 
components are dominant. 


the oscillation wavelength, is the initial phase at the 
orifice and is the decay length. For our liquid jet, 
a is varied in the range of 11 — 15 /rm. The parameters 
Am and Lm from the linear theory are given by following 
formulae: 


Am/Vz = 

27r 

j pa^ 

(2) 

yjmlrn? — 1) ^ 

1 j. , 

Lmhz = 

1 


(3) 

2m(m — 1) V ’ 



where T, p and v are the surface tension, the density and 
the kinematic viscosity of the liquid, respectively, and 
Vz is the speed of the jet. The quantities Amjvz and 
Lral'Vz are the oscillation period and the decay time, re¬ 
spectively. 

Typically the jet is operated at U 2 ~ 10 — 20 m/s. For 
ethanol at 20°C, T = 2.23 x lO-^ N/m, p = 789 kg/m^, 
and V = 1.52 x 10“® m^/s. With these values we obtain 


the period of quadrupolar oscillation Aijvz = 28.0 ^s 
from Eq. ([^. It agrees well with the observed value 
within the experimental error. The decay time ob¬ 
tained from Eq. ^ is = 37.0 ps. Comparing the 

quadrupolar and octapolar modes, we obtain L 2 /L 4 = 6 
and A 2 /A 4 = i/lb- During one period of the quadrupolar 
oscillation, A 2 = 300 — 600/im, the quadrupolar ampli¬ 
tude is halved whereas the octapolar amplitude decays 
by 2 orders of magnitude ® ~ 10“^). Therefore, af¬ 
ter one quadrupolar period, the octapolar mode should 
become hardly noticeable according to the linear theory, 
and thus we expect that the cross-sectional shape of the 
jet at extreme positions Dn with n > 2 [see Fig. 1(c)] 
would be approximated as 

r/a ~ 1-I-? 7 (z) cos 20, (4) 

where ri{z) is identified as the eccentricity e of the cross- 
sectional shape, given by the relation 

^ _ .bjnajor(.^) 

where Lmajoiiz) and Tminor(^) are the lengths of the ma¬ 
jor and minor axes at the observed extreme position, re¬ 
spectively. Note one period of the lowest-order oscillation 
from the nozzle in Fig. 1(c) corresponds to extreme po¬ 
sition D2. 


C. Inconsistency of linear analysis with 
spectroscopic observations 

The above linear analysis cannot account for the spec¬ 
troscopic observations performed on our liquid jet. Be¬ 
fore discussing the details of the inconsistency, let us first 
consider how those observations are made. 

A small z segment of a few micron thickness of the 
liquid column at an extreme position acts as a two- 
dimensional dielectric cavity for the optical wave. The 
boundary of the cavity is given by the cross-sectional 
(perpendicular to the direction of the jet advance) shape 
of the jet. The liquid contains dye molecules. When the 
small z segment of the jet column is excited by a pump 
laser as illustrated in Fig.j^a), the fluorescence from dye 
molecules is enhanced at the cavity resonances as shown 
in Fig. ^b). This enhancement is the cavity quantum 
electrodynamics effect m- In this cavity-modified fluo¬ 
rescence (CMF) spectrum from the micro jet cavity, we 
typically observe 4 — 6 groups of cavity resonances or 
modes. Each mode group is a sequence of resonances 
with a well-defined free spectral range (FSR). 

The CMF spectrum in Fig. [^b) was measured with 
cavity eccentricity e = 0.19 ± 0.02, which was calculated 
by using Eq. ([^ with the measured Tmajor and Lminor 
values from the microscope image. The CMF spectrum 
was taken in a spectral range from 540 nm to 660 nm 
and all cavity modes were identified. The wavelengths of 
those modes were then converted into a normalized size 
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FIG. 3: (a) Experimental setup for CMF spectroscopy (top view). The incident laser is polarized in the direction parallel to 

the jet column, (b) The observed CMF spectrum and identification of five groups (labeled as 1, 2, ..., 5) of cavity modes. The 
spectrum is measured at D4 with a jet-ejection pressure of 1.408 bar. The jet medium is ethanol doped with Rhodamine590 
dye at a concentration of 0.05 mM/L. The mean radius of the jet is a = 14 ± 1 /rm. Positions of the cavity modes are marked by 
vertical bars colored differently for each mode group, (c) Experimentally observed mode evolution diagram when e = 0.19±0.02. 
(d) Mode evolution diagram obtained by the wave calculation for the boundary profile given by Eq. Q with rj = e = 0.19. 
Legends Gl, G2, etc denote cavity mode groups. 


parameter X = 2'Kna/{n^X), where n is the wavelength- 
and temperature-dependent refractive index of the cavity 
medium (ethanol) and no = 1.361, the value of n at 610 
nm and 20° C. Their relative size parameters (labeled as 
X Gap) are then displayed in Fig. [^c) with respect to 
reference size parameters. The reference size parameters 
(displayed on the horizontal axis) are given by an arith¬ 
metic sequence with a mean FSR of all mode groups. The 
diagram obtained in this way as in Fig. [^c) is called the 
‘mode evolution’ diagram [5]. It clearly shows intermode 
interactions, i.e., crossing and avoided crossing of modes. 

These observed modes are then compared with those 
by numerical wave calculation. The wavelength and 
linewidth of a cavity resonance are obtained by solv¬ 
ing the Maxwell equations for a cylindrical dielectric 
medium whose 2-D boundary profile is given by Eq. Q 
(quadrupole deformation only) with 77 = e = 0.19. We 
use the boundary element method for the wave calcula¬ 
tion [in]. 

For visual comparison, the mode evolution diagram ob¬ 
tained from the wave calculation is displayed in Fig. id). 
If Equation Q gives a correct description of the cross¬ 


sectional shape of the jet, the experimental and theoret¬ 
ical mode evolution diagrams should match each other. 
Instead, we find a big discrepancy. In the experiment we 
observe 5 distinct high-Q (quality factor) mode groups 
whereas we only find 4 such groups in the wave calcu¬ 
lation. Moreover, the way the modes interact with each 
other shows a complete discrepancy: the relative spec¬ 
tral positions of all resonance-modes are quite different 
in both cases. 

One may suspect that the uncertainty in the eccen¬ 
tricity might be responsible for such a big discrepancy. 
We find, however, that in order to see 5 high-Q mode 
groups in the wave calculation the eccentricity has to be 
less than 0.14, which is beyond our experimental uncer¬ 
tainty. Moreover, even with such a low eccentricity the 
mode evolution diagram from the wave calculation look 
completely different from the experimental mode evolu¬ 
tion diagram. We then have to conclude that the cross- 
sectional profile of the jet is not given by Eq. (|^. 

What is then the actual cross-sectional shape of the 
jet? The linear theory as an approximate description can 
explain the overall features of the jet such as axis switch- 
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FIG. 4: (a) Top view of our experimental setup and a real 
image of the diffraction pattern, (b) An example of the angle 
width function A0. (c) Reconstructed profile from (b). 


ing and oscillation period. Therefore, it is natural for us 
to assume that the quadrupolar component is truly dom¬ 
inant in the actual boundary profile. There may exist a 
small shape perturbation which might be responsible for 
the different resonance mode interactions. One possible 
candidate for such a small perturbation is an octapolar 
component. The hexapolar component with cos 3d de¬ 
pendence is ruled out by the observed perpendicularity 
of the major and the minor axes or two-fold axial (D 2 ) 
symmetry of the observed cross-sectional shape. 

One may then wonder how such a small perturbation 
can affect the cavity resonances so seriously. Considering 
particle motion in a deformed billiard is helpful to answer 
this question. In quantum chaos, it is a well-established 
fact that there exists a close connection between a clas¬ 
sical trajectory of a particle (in our case, ray trajectory) 
and a wave-mechanical solution (cavity resonances) for a 
given billiard shape [2Q1ET]. It is demonstrated in Ap¬ 
pendix A that a long-time particle or ray dynamics in a 
quadrupolar billiard is very sensitive to a shape perturba¬ 
tion. If experiments are performed with low-Q cavities, 
which is often the case with etched semiconductor mi¬ 
crocavities, long-time dynamics hardly survive to affect 
the cavity resonances much. On the other hand, if we 
deal with high-Q resonances, such as in dynamical tun¬ 
neling and mode evolution, long-time dynamics becomes 
important and thus a small error in the boundary profile 
can lead to a considerable difference [iniiMS]. By this 
reason, the spectroscopy of high-Q resonance modes can 
serve a very sensitive measure of the boundary shape of 
the 2-D optical billiard or the cross-sectional shape of our 
liquid jet. 


III. EXPERIMENTS AND RESULTS 
A. Measurement of the boundary profile 

In our previous studies, we developed an optical tech¬ 
nique to map out the cross-sectional profile of an opaque 
cylindrical object from its forward diffraction patterns 
with a resolution of 0.1% of the mean diameter of the ob¬ 
ject [55]. This technique is particularly useful for fluidic 
columns like a liquid microjet, where the conventional 
electronic imaging techniques such as scanning electron 
microscopy or scanning tunneling microscopy cannot be 
applied. We will use this technique to determine the 
cross-sectional boundary profile of our liquid jet. Below 
we recapitulate its working principle briefly. 

Let us consider the forward diffraction pattern of an 
incident plane wave by an opaque cylindrical object. It 
can be shown that the angular separation A0 between 
two adjacent local minima of the diffracted-wave inten¬ 
sity distribution near the forward direction is given by 
the following simple relation. 

A0 = ^ (6) 

where 2 w is the projected full width of the object in the 
forward direction and A is the wavelength of the incident 
plane wave. We can then reconstruct the cross-sectional 
profile of the object by measuring w at various projection 
angles, covering a range of at least 180° for an object of 
two-fold axial symmetry. Variation of the projection an¬ 
gle is conveniently done by rotating the cylindrical object 
on a rotatable stage for a fixed incident plane wave. The 
details of the boundary profile reconstruction procedure 
is reported in Ref. [55] . 

In our experiment depicted in Fig.|^a), a liquid jet as¬ 
sembly is mounted on a rotary stage and the jet vertically 
ejected is illuminated by a cw laser beam at 532 nm (or 
514 nm in some cases), which is focused by a cylindrical 
lens with a focal length of 7 cm down to a vertical waist 
of 3 ^m on the jet. The diffraction pattern of the laser 
beam by the jet column is then examined on a screen 
located at 2.4 m from the jet. The mean radius of the jet 
column is about 15 ^m while the horizontal waist of the 
laser beam is 2.3 mm, fully covering the jet column hori¬ 
zontally. About 5 — 6 local minima lying within 10° from 
the forward direction are used to obtain a mean value of 
A0 for a given rotation angle 0 of the jet assembly. The 
typical value of the resulting A0 is about 1 degree. 

The requirement of opaqueness is satisfied by doping 
the liquid jet with dye molecules until the optical density 
due to absorption for a laser beam going through the di¬ 
ameter exceeds 6.9. For this, Rhodamine 590 dye, whose 
absorption cross section a at A=532 nm is 3x 10“^^ m^, is 
dissolved in ethanol at a molar density of 23 mM/L. The 
boundary profile measurement is repeated for various ini¬ 
tial conditions such as different ejection pressure, types 
of solvent and mean radii at several critical positions Dn 
(n > 2) on the jet. 
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FIG. 5: B values for various deformations and materi¬ 

als. The black-filled squares and red-filled circles are data 
of ethanol, while the blue- and green-filled triangles are that 
of ethanol-water 50:50 mixture. Numbers are the normal¬ 
ized size parameter X. The averaged B factor (black line) 
is clearly distinguished from the B factor for an ellipse (red 
line a.t B = 0.75) and that for a quadrupole (dotted line at 
B = Q). The fitting error is indicated as a shaded region. 


All boundary profiles we reconstructed exhibit 2- 
fold axial symmetry. The profiles are dominated by a 
quadrupolar component. The second largest component 
is a octapole component. The magnitude of the octapole 
component in almost all cases cannot be explained by 
the linear theory. For example, A0 as a function of 6 in 
Fig .|^b) and its reconstruction profile in Fig.|^c) are ob¬ 
tained at an ejection pressure of 2.6 bar with r = 15 ^m 
and Vz = 16 m/s at D3, which is located at 860 ^m (cor¬ 
responding to about 2 A 2 ) from the orifice. The relative 
magnitudes of the quadrupole and octapole components 
there are 23.9(±0.3)% and 2.6(±0.2)%, respectively. The 
magnitudes of the higher poles are smaller than 0.1%, 
which is within the error range of the present technique. 
If the linear theory were correct, the observed magnitude 
of the octapole component should come from an initial 
relative magnitude exceeding 15000%, which is simply 
impossible. 

Another noticeable point is that there exists a relation 
between the magnitude iji of the quadrupolar compo¬ 
nent and the magnitude 772 of the octapolar component 
as summarized by 

r]2 ~ (7) 


where B is a certain constant to be called ‘i? factor’ be¬ 
low. For an ellipse, from Eq. (A1), we have 5 = 3/4. The 
observed boundary profile of the liquid jet also shows a 
similar relation, but with a different B factor. Figure 
summarizes the B factors extracted from the observed 
boundary profiles under different experimental condi¬ 
tions such as mean radii, types of solvent and eccentric¬ 
ity. Interestingly, they are centered around a common 


value B = 0.42 ± 0.08, if we require 5 to be a con¬ 
stant, despite the different experimental conditions. The 
boundary profile can then be described by the following 
empirical equation: 

r/a ~ 1-b 771 cos20-b 0.42?7 i^ cos40. ( 8 ) 


B. Spectroscopic confirmation of the B factor 

In the preceding section, the spectroscopic observation 
was presented as a counter evidence against the linear 
analysis. Now, for the boundary profile given by the em¬ 
pirical formula Eq. (|^ with 771 = 0.185±0.003, the exper¬ 
imental mode evolution diagram in Fig. j^c) is in good 
agreement with the wave calculation result. This 771 value 
corresponds to an eccentricity e of 0.182 ±0.003, which is 
consistent with the measured eccentricity of 0.19 ± 0.01 
for our liquid jet within the experimental error. Besides, 
in Fig. five more mode diagrams are presented for a 
wide range of 771 . 

One can evaluate the tolerance in the B factor value 
of 0.42 by examining the difference between the observed 
mode evolution diagram and the calculated one as a func¬ 
tion of B around the empirical value of 0.42. For quan¬ 
titative analysis, we pay attention to the avoided cross¬ 
ing gap between mode groups 2 and 4, which are de¬ 
noted by G2 and G4, respectively, in Figs. and The 
G2-G4 avoided crossing gap for 771 = 0.185 is calculated 
from the wave equation for various B factor values from 
0.37 to 0.45 in Fig. [^a). We find that the experimen¬ 
tal avoided crossing gap is well fit by the B factor value 
of 0.414 ± 0.011. This is also supported in Fig. [^b), 
where both wave calculations with B=0.41 and 0.42 well 
reproduce the experimental mode evolution diagram for 
771=0.185. 

We can now regard the empirical formula Eq. § as 
a correct description of the actual cross-sectional bound¬ 
ary profile of our liquid jet. The remaining question is 
then how we can derive the empirical formula from the 
hydrodynamics theory m by going beyond the linear 
approximation. 


IV. NONLINEAR ANALYSIS 
A. Bohr’s analysis 

In 1909, Niels Bohr at the age of 23, before inventing 
the hydrogen model, considered the liquid jet problem 
including nonlinearity m- Conceiving the limitation of 
the linear approximation in the previous study by P. O. 
Pederson [28], he derived an expression for the surface 
oscillation of a non-viscous liquid jet with a few lowest- 
order components while taking the nonlinear interactions 
among them into account |29j . 

When we rewrite his result in terms of our notations, 
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FIG. 6: Matching the mode evolution diagrams between experiment (denoted by ‘Exp’) and wave calculations (denoted by 
‘Cal’). Experimental conditions are expressed as ‘(jet pressure in bar, extreme position Dn, jet temperature in ‘’C)’. Jet mean 
radius a for each case is determined as a fitting parameter, (a) Experiment at (1.250, D5, 16.9) and the wave calculation with 
rji — 0.095. a=13.79±0.01 fim. (b) Experiment at (1.402, D5, 16.9) and the wave calculation with rji = 0.105. a=13.71±0.01 
fim. (c) Experiment at (1.326, D5, 21.4) and the wave calculation with r/i = 0.130. a=13.64±0.01 fim. (d) Experiment at 
(1.188, D4, 20.3) and the wave calculation with rji = 0.160. a=13.61±0.01 /rm. (e) Experiment at (1.408, D4, 19.5) and 
the wave calculation with rji — 0.185. a=13.55±0.01 fim. (f) Experiment at (1.536, D4, 19.1) and the wave calculation with 
rji = 0.200. a=13.50±0.01 /rm. The data in (e) is the same as the data in Fig. 3(c). 


we obtain 


r/ao 


1 1 2 1 2 
1 — -ry — -r] 


+ -77^ cos 40 + 
4 


cos 2 kz + rj cos 20 cos kz 

-77^ cos 40 cos2fcz, 

6 ' 


where k is the wave vector in the propagation direction, 
and ao is the mean radius when rj = 0. Considering only 
nodal {cos kz = cos2kz = 1) and anti-nodal {cos kz = 
— l,cos2kz = 1) positions, and taking the major axis of 
the shape as the 0 axis in the r-0 coordinates, Eq. ^ is 


( 9 ) 
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FIG. 7: (a) G2 — G4 avoided crossing gaps from the experiment (open circles with error bars) is compared with those from 

the wave calculation with the B factor varied from 0.37 to 0.45 forn 771 =0.160, 0.185, 0.200. (b) Mode evolution diagrams for 
B = 0.41 and B = 0.42 for t?i = 0.19. 


simplified as 

1 5 

r/ao = 1 —+ rj cos 29 +cos 46. (10) 

From the flux conservation of the jet with its negligible 
speed change, we find 

a/oo~l-^. (11) 

Including this effect, the empirical formula Eq. (|^ can 
be written as 

r/ao ~ 1 — + 77 cos20 + 0.42?7^ cos40, (12) 


If we limit ourselves to the region where the surface 
oscillations have not decayed sufficiently, i.e., at several 
oscillation wavelengths away from the orifice, the change 
of velocity due to gravity is estimated to be negligibly 
small for our jet (less than 10 “^ of the initial velocity). 
The velocity change due to the flux conservation, how¬ 
ever, amounts to 1 % of the initial velocity, in the same 
order of magnitude as the octapole component in the 
typical boundary profile. This is not small enough to be 
neglected, and it may affect the boundary profile. Conse¬ 
quently, we cannot rule out the possibility that the agree¬ 
ment between the experimental result Eq. ([^ or Eq. (12) 
and Bohr’s analysis Eq. (10) might be accidental under 
the particular conditions that the experiment was done. 


which is almost the same as Bohr’s result. Only differ¬ 
ence is the numerical factor 0.42 compared to a rational 
number 5/12 = 0.416. The calculation of Bohr’s gives 
the theoretical B factor of 5/12, which agrees well with 
our empirically B factor within the experimental error. 


In order to resolve this problem, we have extended 
Bohr’s inviscid 2-D analysis into 3-D by including the 
effects of gravity and flux conservation. We have then 
obtained the coefficient attached to rj^ cos 49 as 


B. Limitation of Bohr’s analysis and our extension 


In Bohr’s analysis, the ^-dependence is treated quite 
simply: the surface motion is regarded as a two- 
dimensional oscillation and its time dependence is sim¬ 
ply replaced by the z-dependence with an assumption 
that the fluid velocity is constant. This treatment is by 
no means rigorous by the following reasons. First, un¬ 
der gravity, the fluid velocity is not constant. Particu¬ 
larly, if the jet is ejected vertically upward like in our 
experimental setup, the velocity slows down. Second, 
the conservation of the flux dictates the velocity has to 
change periodically. The sectional area given by Eq. (10) 
is 7r[l -I- \rf{4 — cos2kz)], having z dependence. There¬ 
fore, the velocity also has to change as z in order keep 
the flux conserved. 


B =^- ^ 0 4^7 _ o. 07 ao"fc 2 . (13) 


The detailed derivation is presented in Appendix B. The 
term containing (oo/c)^ is very small in our experiment: 
O.OToo^fc^ ~ 0.42T/ paoVz^ ~ 0.8 [w?/s^\/{v ~ 10 “^ 
considering that the velocity Vz is typically 10 — 20 m/s. 
Experimentally, such a small correction is obscured by 
our experimental error. Therefore, our calculation con¬ 
firms that the nonlinear analysis by N. Bohr is adequate 
for describing the motion of our liquid microjet. How¬ 
ever, if oofc is comparable to or larger than unity such as 
in the case of garden hose or home faucet, the adequate¬ 
ness of Bohr’s analysis is not guaranteed and we expect 
that our extension should be applied. 
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V. CONCLUSION 

We have measured the cross-sectional shapes of a sur¬ 
face oscillating liquid jet at its extreme positions by us¬ 
ing the optical forward diffraction method and deter¬ 
mined the magnitude of their octapole components. The 
normalized magnitude of the octapole component was 
found to be empirically given by Brj^ with rj the magni¬ 
tude of a quadrupole component and B = (0.42 ± 0.08). 
We also measured the resonance mode spectra of a two- 
dimensional cavity formed by a cross-sectional segment 
of the liquid jet at various extreme positions. The ob¬ 
served spectra agreed well with wave calculation assum¬ 
ing a proportionality coefficient B = 0.415±0.010. Com¬ 
bining these two results we were able to confirm the long¬ 
standing prediction of Niels Bohr’s in his nonlinear the¬ 
ory of fluidic surface oscillation for the first time. 

Our finding can be utilized in precise measurement of 
surface tension especially when the material is in a non¬ 
equilibrium state such as being sprayed, spun, or blown. 
Moreover, since our results provide an accurate shape 
formula for the surface oscillation of a fluidic jet, one can 
match experimental data such as mode spectrum and 
output directionality from a given fluidic microjet cav¬ 
ity more accurately with theoretical ray and wave sim¬ 
ulations. One can also predict mode dynamics for arbi¬ 
trary cavity deformation in a fluidic microjet cavity in 
advance, and thus can help designing various photonics 
experiments based on intermode interactions. One exam¬ 
ple is experimental verification of an enhanced gain due 
to the diverging Petermann factor near an exceptional 
point arising from intermode interactions. Fluidic micro¬ 
jet cavities can be incorporated into optofluidic systems 
[4] as optical resonators whose mode characteristics can 
be accurately designed with our results. Such systems 
can be a versatile platform for biosensors [30] j optical 
filters and mode switchers |3T] . 

This work was supported by the Korea Research Foun¬ 
dation (Grant No. 2013R1A1A2010821). 


Appendix A: Effect of a small perturbation in RAY 
dynamics 

It is known that non-regular boundary shapes such as 
quadrupole can induce chaotic dynamics in the classi¬ 
cal particle description and also result in the interactions 
among resonance modes in the wave description. In ad¬ 
dition, it is known that even quite a small perturbation 
can change the particle dynamics significantly in a par¬ 
tially chaotic system. Considering the close connection 
between classical dynamics and wave solutions, we ex¬ 
pect a small perturbation of a right type may explain the 
above discrepancy with the spectroscopic observations. 

A regular billiard which has a profile similar to a 
quadrupole is an elliptical billiard, whose boundary can 
be expressed as an infinite series of multipolar compo- 



FIG. 8: Comparison of an elliptical billiard (black curve), a 
semi-elliptical billiard (blue curve) and a quadrupolar billiard 
(red curve) with the same eccentricity e = 0.16 and the same 
enclosed area. When a particle is reflected off the boundary, 
its angular position is denoted by 0 with an reflection (or 
incident) angle x- 

nents as 

3 5 

rja^ I ± 77 cos 20 ± - 77 ^ cos 40 ± - 77 ^ cos 60 • • • (Al) 

4 8 

The shape difference between a pure quadrupole and the 
ellipse with the same eccentricity is quite small, amount¬ 
ing to at most 1/100 of the mean diameter 2 a in the de¬ 
formation range of 0 < 77 < 0.2 as seen in Fig.[^ Beyond 
77 = 0 . 2 , the convexity of quadrupole is broken like peanut 
due to the contraction along the minor axis. However, 
internal ray dynamics in the elliptical billiard are quite 
different from those of the pure quadrupolar billiard. It 
is because a small initial difference in particle trajectories 
can grow exponentially later on in the presence of chaos. 

Classical particle dynamics is conveniently analyzed in 
Poincare surface of sections (PSOS), which is commonly 
used to visualize the particle motion in a 2 -dimensional 
billiard. PSOS is a phase-space diagram which is drawn 
iteratively in terms of reflection position ( 0 ) and reflec¬ 
tion angle (x) of a particle upon bouncing off the bound¬ 
ary. Regular or quasi-periodic orbits are represented by 
concentration of dots around localized substructures such 
as islands or curves. Chaotic motion has no such con¬ 
centration of dots, rather randomly but uniformly filling 
some regions called chaotic sea in the phase space. 

Figurej^a) shows the PSOS of a pure quadrupole while 
Fig. I^b) shows that of an ellipse, both drawn in terms 
of 0 and siny for e = 0.16. Almost the whole PSOS 
in Fig. ia) is chaotic while the PSOS in Fig. [^b) con¬ 
sists of many periodic and quasi periodic orbits without 
chaos. Since it is practically impossible to fabricate a bil- 
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’to 



FIG. 9: Poincare surface of sections for various billiards, (a) A pure quadrupole given by 1 + 0.16 cos 20, (b) an ellipse and 
(c) a semi-ellipse, both with an eccentricity e = 0.16. 


hard with a perfectly elliptical or perfectly quadrupolar 
boundary, it is also meaningful to consider the PSOS of a 
truncated semi-ellipse with the same eccentricity for com¬ 
parison as in Fig.j^c). The profile of a semi-elliptical bil¬ 
liard is obtained by truncating the expansion in Eq. (A1) 
beyond the second order of rj. Though chaotic sea is 
present in considerable regions of the phase space, we 
still have curve-like or island-like structures correspond¬ 
ing to regular trajectories. From these considerations, we 
find that even quite a small difference in the boundary 
profile can change the particle dynamics significantly. In 
particular, degradation of the regular dynamics is strik¬ 
ing when an ellipse is replaced with a semi-ellipse with 
the same eccentricity. 


Appendix B: Improved nonlinear analysis 

In order to derive the expression for the B factor more 
rigorously, we have modified the nonlinear calculation of 
Niels Bohr to include the z-dependence up to the second 
order of rj. We include the viscous decay in terms of an 
exponential factor multiplied to an inviscid surface os¬ 
cillation mode. Then the surface motion can be treated 
by introducing the potential (p for the velocity v, satis¬ 
fying V(/> = V. Master equations of motion derived from 
the Bernoulli’s principle under the stationary condition 
are then written as follows in the cylindrical coordinates 
(r,0,z): 


V"<^ = 0, 


(Bl) 


pa\Ri R 2 J 2\drj 2r^\d0j 2\dzj 

(B2) 


d(j) 1 d(j) dr 


d(j) dr 


dz dz 


= 0 , 


(B3) 


where C is a constant and Ri and R 2 represent two 
principal radii of curvature of the jet surface. Equa- 


the Bernou 


tions (Bl )-(B3) account for the irrationality of the fluid. 


li’s principle and the condition for unbroken 


surface, respectively. The sum of the reciprocals of the 


principal radii is just the trace of the shape operator in 
differential geometry: 


gp gp 
i?i i?2 


Q d'^r 

-r^ -h 

dz‘^ 

d'^r ( dr 


2 r 


902 V dz 


dr dr d^r 
do dz dOdz 
dr 


■ 


-rHp) +r= 

. dz, 


^/dry d'^r 


(B4) 


The surface motion is decomposed by the orders of mag¬ 
nitude in rj: 


r/ao = -b 
(j) = -b rjy'^ + 


(B5) 

(B6) 


Then Eq. B4 becomes 


Op , gp , 

Ri^ R 2 ~ 

1 , 


-bgp 


96*2 


az2 


4iz(i) 


( 2 ) d'^F^'^'> ( 




+gp 


( 96<2 


do 


^( 1 ) 

dz 


dz"^ 


(B7) 


The zeroth motion corresponds to an upward cylindrical 
flow of a circular cross section. The gravitation causes 
the vertical velocity to decrease as the height: = v^'^ — 

2 gz, where g is the gravitational acceleration and Vz is 
the initial velocity. At the same time the radius has to 
increase to conserve the flux. Therefore, and 
are given by 


= 1 -b 


gz 

2 u 2 


(B8) 
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<^( 0 ) = 


gr 


I ^ 2 

2 vz 


(B9) 


The r-dependence of the velocity potential is required 
to satisfy the Laplace equation. With these expres¬ 
sions substituted in Eq. (B2), we can identify C = 
T/{pao) + u^/2 — gz. However, the terms containing the 
gravitational acceleration g are too small to make notice¬ 
able contribution. They will be thus neglected from now 
on. 

The motion corresponding to the first-order terms is 
quite similar to the motion from the linear analysis. The 
master equations are: 


= 0 , 


T 

pao 


p(A) 


a2j^(i) 

502 


ao 


,^ 2 ^( 1 ) 

5z2 


dz 


d(j)W 

a^Vz ^:.— = 0. 


dr 


dz 


(BIO) 

= 0 , 

(Bll) 

(B12) 


A solution of the Laplace equation in the cylindrical co¬ 
ordinates can be written as a linear sum of Bessel’s func¬ 
tions J„. Therefore above equations are solved by sub¬ 
stituting 


= cosnOcoakz 

(j)^d — MnJn{ikr) cos nO sink z 


(B13) 


where n = 1 , 2 ,3,... and Mn is a coefficient to be deter¬ 


mined. Equations (Bill and (B12) then yield 


M„ = - 


pakvzJniikao) 


{v? + k^ao^ — 1) = 


laoVz 


J'niikao)' 

(B14) 

The factor (kao) in the Bessel’s function is noticeable. It 
is interpreted as an index number representing the ratio 
of the mean radius to the z-axial wavelength. A circular 
cylinder is t he cas e when fcoo —>■ 0 . 

Equation (B14) is further simplified when fcog <C 1. 


Using the approximated expression for the Bessel’s func¬ 
tion, the above relation can be approximated as 


(fcoo)^ = 


pao^Vz 


n(n^ — 1) 


3n — 1 ,, , 2 

1 + - —{kao) 


2 n{v? — 1) 


or by solving for fcoo we get 


(B15) 


ka ~ yJn{rF — 1) 


T 


1/2 r 


1 -f 


3n - 1 / T 


paov^j [ 4 \paov'^, 

(B16) 

When the term proportional to (kao)^ is negligible, the 
equation is reduced to the wave number formula derived 
in the linear analysis, Eq. ([^. When n = 1 (dipolar mo¬ 
tion), Eq. (B16) is satisfied if and only if fc = 0. The 
dipolar motion corresponds to a constant translation of 
the center of cross-sectional mass. It can be neglected 


as far as surface oscillation is concerned. Therefore, the 
quadrupolar (n = 2 ) motion is the lowest oscillatory mo¬ 
tion. 

Besides the Laplace equation, the second ordered equa¬ 
tions are more complicated. 


T (f<-, + IPF 

pao V 502 


,52J7(2)\ 


-f 


T 

pao 


-2 




dz 

Q 2 fW 




-ao 


\ dz J 


\ de J 5z2 

\ 5r y ao2 \ 50 / 




■ AvzF^^'f 


d(j)'d'> 

dz 


■ 2 aoVzF(^'> 


52 (^( 1 ) 
dzdr 


= 0 , 


(B17) 


d4>^‘^^ 

dr 


aoVz 


dF^^'i 2 5y>(i) 5 F(i) 


dz 


dr a do dO 
5^(1) dF^d 


- AaoVzF^^'>^ -^“ 0 ^-^ 

dz dz dz 


(B18) 


By substituting the expr ession s in Eqs. ( B13) and (B16) 
into F^^'> and Eqs. (B17l and (B18) become 


- ( -h 


52 ^( 2 ) 2 52 ^( 2 ) \ 


{n^ — 1) 2(n2 — l)ao^fc^ 

T, 2 112 


cos^nOcos^kz 


T 

pao 

pa^z t/y-r„2gjj^2^^^Qg2^^ 

1 dz 

— (2 — 4ao^/c^)cos^n0cos^fcz — ao^fc^cos^n0sin^fcz] 




+ ^^ao'^k^Vz^cos^nOsin^kz + 


T 


pao 


n?{ri^ — 1) 

ao2fc2u^2 


-f 


2rr{'nr — 1 ) rrao^k 


21,20 


sin^n 0 sin^fcz 


Vz" Vz 

— A{v? - 1 - 1 - ao^k'^)cos^n9cos^kz 
2pao 


T 


-ao 


^k^Vz'^cos^nOcos^kz = 0, 


(B19) 




- aoVz- 


5F(2) 

dz 


^(n^ - 1 - 1 - ao^k^)src?n9 sin 2 kz 


dr 

T 1 
pao aokvz 

-n^(n^ — 1 -I- Oo^A:^)cos^n0sin2A:z = 0. 

(B20) 

The terms comparable to or smaller than ao'^k'^T/[pao) 
have been neglected since the sum of them would con¬ 
tribute the final result by no more than 0 . 01 %. 


T aok 2/ 2 1 , „ 2j2\ _ 


pao Vz 
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Now we are ready to evaluate the second ordered mo¬ 
tion . The velocity p otent ial must satisfies both 
Laplace equation and Eq. (B20|. Therefore, and 
can be written as 


77^ cos40, is then obtained from 6 ii{z) with n = 2 (oc- 
tapole component) and cos 2kz = 1 (nodal and anti-nodal 
positions): 


= Jo{i2kr)joiz) + J2nii‘2kr)cos2n9j2niz), 

= (5 o( 2:)-f cos 2n0(52n(-z)- (B21) 


Bn=2 — 2 '^4(2)|cos2fc2=l 


5 

12 


232 

3240 


2 7 2 

ao k 


~ 0.417-O.OToo^fc^. 


(B23) 


The function Sq{z) = sin^ kz is related to the area con¬ 
servation. Carrying on the calculation further yields the 
expression for S 2 niz). 


which is Eq. (13) quoted in Sec. IV B 


2n^ + 5n^ — 2n — 2 n'^ + n + 1 272 

= -4(4n2 - 1)-+ 4n(4n2-l)“° ^ 

—2n^ -I- 7n? +2n — A 

^ [ 4(2n2 -p 1) 

—40n^ -|- 22n'^ -|- 23n^ -|- 17n^ -I- lln -|- 1 ,,9 

-H - 9 - ao k^ 

An{2n -I- l){2n‘^ -|- 1 )^ 

X cos2fcz. (B22) 


The expression for Bn= 2 , the coefficient attached to 


[1] D. Psaltis, S. R. Quake, and C. Yang, Developing 
Optofluidic Technology through the Fusion of Microflu- 
idics and Optics, Nature 442, p. 381 (2006). 

[2] C. Monat, P. Domachuk, and B. J. Eggleton, Integrated 
Optofluidics: A New River of Light, Nat. Photonics 1, p. 
106 (2007). 

[3] V. R. Horowitz, D. D. Awschalom, and S. Pennathur, 
Optofluidics: Field or Technique?, Lab Chip 8, p. 1856 
(2008). 

[4] Y. Fainman, L. P. Lee, D. Psaltis, and C. Yang, Optoflu- 
tdics (McGrawHill, New York, 2010). 

[5] J. Yang, S. Moon, S.-B. Lee and J.-H. Lee, K. An, De¬ 
velopment of a Deformation-Tunable Quadrupolar Micro- 
cavity, Rev. Sci. Instrum. 77, 083103 (2006). 

[6] J. U. Nockel, A. D. Stone, G. Chen, H. L. Grossman, 
and R. K. Chang, Directional Emission from Asymmetric 
Resonant Cavities, Opt. Lett. 21, p. 1609 (1996). 

[7] S.-B. Lee, J. Yang, S. Moon, J.-H. Lee, and K. An, J.- 
B. Shim, H.-W. Lee, S. W. Kim, Universal Output Di¬ 
rectionality of Single Modes in a Deformed Microcavity, 
Phys. Rev. A 75, 011802(R) (2007). 

[8] S.-B. Lee, J. Yang, S. Moon, J.-H. Lee, and K. 
An, Chaos-Assisted Nonresonant Optical Pumping of 
Quadrupole-Deformed Microlasers, Appl. Phys. Lett. 90, 
041106 (2007). 

[9] S.-B. Lee, J. Yang, S. Moon, S.-Y. Lee, J.-B. Shim, S. W. 
Kim, J.-H. Lee, and K. An, Quasieigenstate Evolution 
in Open Chaotic Billiards, Phys. Rev. A 80, 011802(R) 
(2009). 

[10] J. Yang, S.-B. Lee, S. Moon, S.-Y. Lee, S.-W. Kim, T. T. 
A. Dao, J.-H. Lee, and K. An, Pump-Induced Dynamical 
Tunneling in a Deformed Microcavity Laser, Phys. Rev. 
Lett. 104, 243601 (2010). 

[11] S. E. Bechtel, M. G. Forest, N. T. Youssef, H. Zhou, The 


Effect of Dynamic Surface Tension on the Oscillation of 
Slender Elliptical Newtonian Jets, J. Appl. Mech. 65, p. 
694 (1998). 

[12] D. Narita, M. Baba, and K. Ohtani, Three-Dimensional 
Shape Measurement of a Transparent Object Using a 
Rangefinding Approach, Instrumentation and Measure¬ 
ment Technology Gonference, IMTC 2003. Proceedings 
of the 20th IEEE 2, p. 1022 (2003). 

[13] Lord F. R. S. Rayleigh, On the Instability of Jets, Proc. 
Lond. Math. Soc. 10, p. 4 (1878). 

[14] Lord F. R. S. Rayleigh, On the Capillary Phenomena of 
Jets , Roy. Soc. Proc. 29, p. 71 (1879). 

[15] Sir H. Lamb, Hydrodynamics (Dover Publications, New 
York, 1945). 

[16] N. Bohr, Determination of the Surface-Tension of Water 
by the Method of Jet Vibration, Philos. Trans. Roy. Soc. 
Lond. 209, p. 281-317 (1909). 

[17] T. V. Kasyap, D. Sivakumar, B. N. Raghunandan, Flow 
and Breakup Characteristics of Elliptical Liquid Jets, Int. 
J. Multiphase Flow 35, p. 8 (2009). 

[18] P. R. Berman, Cavity Quantum Electrodynamics (Aca¬ 
demic Press, Boston, 1994). 

[19] J. Wiersig, Boundary Element Method for Resonances in 
Dielectric Microcavities, J. Opt. A: Pure Appl. Opt. 5, 
53 (2003). 

[20] J. U. Nockel and A. D. Stone, Ray and Wave Chaos in 
Asymmetric Resonant Cavities, Nature 385, 45 (1997). 

[21] J.-B. Shim, S.-B. Lee, S. W. Kim, S.-Y. Lee, J. Yang, S. 
Moon, J.-H. Lee, and K. An, Uncertainty-Limited Turn¬ 
stile Transport in Deformed Microcavities, Phys. Rev. 
Lett. 100, 174102 (2008). 

[22] S. G. L. Harald, N. B. Rex, H. E. Tureci, R. K. Ghang, 
A. D. Stone, T. Ben-Messaoud, and J. Zyss, Dramatic 
Shape Sensitivity of Directional Emission Patterns from 













13 


Similarly Deformed Cylindrical Polymer Lasers, J. Opt. 
Soc. Am. B 21, p. 923 (2004). 

[23] T. Tanaka, M. Hentschel, T. Fukushima, and T. 
Harayama, Classical Phase Space Revealed by Coherent 
Light, Phys. Rev. Lett. 98, 033902(2007). 

[24] S.-B. Lee, J. Yang, S. Moon, S.-Y. Lee, J.-B. Shim, S. W. 
Kim, J.-H. Lee, and K. An, Observation of an Exceptional 
Point in a Chaotic Optical Microcavity, Phys. Rev. Lett. 
103, 134101 (2009). 

[25] S. Shinohara, T. Harayama, T. Fukushima, M. Hentschel, 
T. Sasaki, and E. E. Narimanov, Chaos-Assisted Di¬ 
rectional Light Emission from Microcavity Lasers, Phys. 
Rev. Lett. 104, 163902 (2010). 

[26] S. Moon, J. Yang, S.-B. Lee, J.-B. Shim, S. W. Kim, J.- 
H. Lee, and K. An, Nondestructive High-Resolution Soft- 
Boundary Profiling Based on Forward Shadow Diffrac¬ 
tion, Optics Express 16, 11007 (2008). 

[27] It should be pointed out that the hydrodynamic master 
equations treated with Poincare-Lindstedt method fail 


to include all of the nonlinear terms required in Bohr’s 
nonlinear theory. For example, see J. A. Tsamopoulos 
and R. A. Brown, Non-Linear Oscillations of Inviscid 
Drops and Bubbles, J. F. Mech. 127, p. 519 (1983). 

[28] P. O. Pederson, On the Surface-Tension of Liquids In¬ 
vestigated by the Method of Jet Vibration, Philos. Trans. 
Roy. Soc. Lond. A 207, p. 341 (1908). 

[29] Since Bohr’s intention was to measure the surface ten¬ 
sion, no attention was paid to the boundary profile at 
extreme positions by himself. Likewise, the aspect of his 
work on the cross-sectional boundary profile has been 
completely unrecognized. 

[30] A. M. Armani, R. P. Kulkarni, S. E. Fraser, R. C. Flagen, 
and K. J. Vahala, Label-Free, Single-Molecule Detection 
with Optical Microcavities, Science 317, p. 783 (2007). 

[31] V. S. Ilchenko and A. B. Matsko, Optical Resonators with 
Whispering-Gallery Modes-Part II: Applications, IEEE 
J. Sel. Top. Quantum Electron. 12, p. 15 (2006). 



